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SUMMARY 

Films  of  magnesium  fluoride  are  found  to  be  under  large 
tensile  stress  following  exposure  to  water  vapour.  It  is 
shown  that  this  stress,  which  varies  reversibly  in  magnitude 
with  the  adsorption  equilibrium,  can  be  attributed  to  electro¬ 
static  dipole  interaction  between  the  adsorbed  water  molecules. 

The  effect  is  not  confined  to  magnesium  fluoride,  but  can 
occur  in  other  materials  of  an  appropriate  porous  structure. 

Its  relevance  to  the  stability  of  optical  coatings  is 
discussed,  together  with  some  other  applications. 
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1.  INTRODUCTION 

In  the  course  of  an  investigation  into  the  stress  behaviour  of  vacuum 
evaporated  films  it  was  observed  that  when  freshly  deposited  layers  of  magnesium 
fluoride  were  exposed  to  the  atmosphere,  they  immediately  showed  evidence  of 
very  substantial  tensile  stress.  The  phenomenon  is  not  confined  to  MgF2  films 
only,  and  in  view  of  this  more  general  importance,  particularly  in  relation  to 
the  stability  of  optical  coatings,  it  was  decided  to  examine  the  effect  in  more 
detail.  The  purpose  of  this  paper  is  to  describe  the  experimental  observations, 
and  to  propose  a  model  for  the  mechanism  leading  to  the  production  of  the 
tensile  stress. 


2.  EXPERIMENTAL  ARRANGEMENT  AND  OBSERVATIONS 


The  experiments  were  made  on  MgF2  films  deposited  onto  thin  mica  substrates 
of  4  mm  width,  clamped  at  one  end  so  as  to  form  cantilevers  of  20  mm  free 
length,  and  stress  measurements  were  carried  out  by  the  well  known  method  of 
optically  observing  the  deflection  5  of  the  free  end  of  the  substrate(ref .1) . 
With  the  usual  simplifying  assumptions  inherent  in  this  technique,  and  neglect¬ 
ing  the  Poisson  ratio  correction,  the  film  stress  per  unit  width  is  given  by 
(ref. 2) : 


4ED2  5 
3L2Tf 


where  E  = 


D  = 
L  = 


Young's  Modulus  of  the  film 
thickness  of  the  substrate 
length  of  cantilever 
film  thickness 
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Cleaved  mica  strips  with  thicknesses  in  the  range  from  0.016  mm  to  0.025  mm 
were  used  as  substrates.  For  each  of  these  strips  E  was  determined  by  observ¬ 
ing  the  beam  deflections  produced  by  known  weights,  the  measurements  yielding 
an  average  value  of  1.8  x  1010  kg/m  ,  with  little  variation  from  sample  to 
sample. 

The  film  depositions  and  stress  measurements  were  performed  in  the  apparatus 
shown  schematically  in  figure  1.  A  liquid  nitrogen  trapped  oil  diffusion  pump 
permitted  the  experimental  chamber  to  be  evacuated  to  a  base  pressure  of  about 
5  x  10-6  Torr,  following  which  the  system  could  be  isolated  from  the  pump  by  a 
baffle  valve.  Gas  could  then  be  admitted  to  any  desired  pressure  p,  measured 
on  a  membrane  manometer.  During  this  phase  of  the  procedure  condensible 
components  of  the  filling  gas,  such  as  water  vapour,  could  be  removed  from  the 
experimental  volume  by  means  of  a  large,  liquid  nitrogen  cooled  freezing  pocket. 

The  deflection  of  the  cantilever  was  observed  through  a  window,  using  an 
optical  filar  micrometer,  capable  of  detecting  movement  down  to  about  0.01  mm. 

For  reasons  of  experimental  convenience  the  MgF2  films  of  about  4000  A  thick¬ 
ness  were  deposited  with  the  vapour  stream  incident  at  an  angle  of  roughly  45° 
to  the  substrate,  which  was  at  ambient  temperature. 

Earlier  findings  on  MgF2  films (ref . 3,4)  produced  on  glass  or  quartz  sub¬ 
strates,  have  usually  shown  the  films  to  possess  an  intrinsic  tensile  stress. 

In  contrast  to  this,  the  present  films,  laid  down  on  mica,  were  invariably  under 
a  compressive  stress,  typically  of  the  order  of  109  kg/m2 .  The  corresponding 
beam  deflection  remained  quite  stable  as  long  as  the  film  was  kept  in  a  good 
vacuum,  but  immediately  after  exposure  to  the  atmosphere  the  strip  began  to  move 
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in  the  direction  opposite  to  its  initial  deflection.  This  response  was  rapid, 
and  under  suitable  circumstances  it  continued  well  beyond  the  original  "neutral" 
position,  adopted  by  the  strip  prior  to  deposition  of  the  film,  thus  indicating 
that  a  net  tensile  stress  had  been  generated  by  exposure  to  the  atmosphere. 

The  effect  was  reversible  in  the  sense  that  when  the  experimental  chamber  was 
re-evacuated,  the  beam  again  took  up  the  position  corresponding  to  the  intrinsic 
compressive  stress  found  immediately  after  deposition,  and  moreover  the  sequence 
of  air  admission  and  evacuation  could  be  cyclically  repeated  with  good 
reproducibility.  We  conclude  therefore  that  the  original  compressive  stress 
remained  preserved  at  a  constant  level  throughout,  but  that  superimposed  on  it 
was  a  variable  tensile  component,  that,  as  we  shall  see  below,  depended  in 
magnitude  on  the  properties  of  the  ambient  atmosphere. 

A  series  of  experiments  established  unequivocally  that  admission  of  dry 
permanent  gases  such  as  oxygen,  nitrogen,  hydrogen  or  argon  was,  at  pressures  up 
to  one  atmosphere,  quite  ineffective  in  producing  the  tensile  component  of 
stress,  but  that  the  presence  of  water  vapour,  even  at  pressures  well  below  one 
Torr,  was  sufficient  to  cause  a  readily  observable  tensile  force,  indicating  that 
this  was  in  fact  the  only  atmospheric  constituent  to  contribute  significantly  to 
the  effect.  Accordingly  the  subsequent  work  was  focussed  on  determining  the 
role  played  by  H2O.  To  this  end  moist  laboratory  air  was  introduced  into  the 
chamber  at  pressures  of  several  hundred  Torr,  and,  with  the  baffle  valve  closed, 
its  water  vapour  component  was  condensed  in  the  freezing  pocket.  After  the 
system  had  reached  stability,  the  liquid  nitrogen  was  removed  from  the  pocket, 
and  the  deflection  of  the  beam  was  observed  as  a  function  of  the  changing 
pressure  as  the  condensed  water  vapour  was  released  into  the  chamber.  Separate 
experiments  confirmed  that  the  beam  position  was  not  affected  by  stray  thermal 
effects  due  to  the  presence  or  absence  of  the  coolant. 

The  response  of  the  system  was  greatly  accelerated  if,  after  condensation  of 
the  water  vapour,  the  baffle  valve  was  opened  for  a  period  long  enough  to 
evacuate  the  non-condensible  gases  from  near  atmospheric  pressure  to  a  level 
below  1  x  10"5  Torr,  thus  facilitating  the  subsequent  diffusion  of  the  vapour 
through  the  apparatus.  This  procedure,  which  left  H20  as  the  only  significant 
constituent  of  the  residual  atmosphere,  was  considered  the  preferred  experimental 
method,  although  the  ultimate  result  of  the  observations  was  not  affected  by  the 
presence  of  permanent  gases  at  high  partial  pressure. 

Observations  on  a  number  of  films  confirmed  that  the  stress  behaviour  was  very 
similar  in  all  cases.  As  an  example  figures  2(a)  apd  2(b)  summarise  two 
typical  experimental  runs,  both  involving  the  same  film,  but  with  a  different  air 
filling  used  for  each  measuring  sequence.  The  curves  represent  the  beam 
deflection  8  and  the  water  vapour  pressure  p  as  functions  of  the  time  t,  measured 
from  the  instant  of  removal  of  the  coolant  from  the  freezing  pocket.  Because  of 
small,  uncontrolled  differences  in  the  detailed  manner  in  which  liquid  nitrogen 
was  applied  and  withdrawn,  the  pressure-time  dependence  varied  somewhat  between 
runs.  In  particular  figure  2(b)  gives  an  example  of  a  pressure  overshoot  that 
was  occasionally  noticed  before  dynamic  adsorption-desorption  equilibrium  was 
established  in  the  experimental  chamber.  It  is  seen  that  the  resulting 
pressure  maximum  is  associated  with  a  corresponding  maximum  in  the  8  curve, 
demonstrating  not  only  the  reversible  nature  of  the  tensile  stress,  but  also 
indicating  roughly  the  speed  with  which  the  stress  responds  to  pressure  changes. 

Although,  as  we  have  seen,  the  variations  of  8  and  p  with  time  are  different 
in  figures  2(a)  and  2(b),  a  plot  of  5  versus  p,  using  data  from  both  runs,  is 

the  smooth  curve  of  figure  2(c),  proving  the  existence  of  a  definite  5  -  p 

relation. 

Since  on  the  one  hand  H2O  appears  to  be  the  only  atmospheric  constituent 
involved,  and  since  on  the  other  hand  the  water  molecule  is  also  the  only  one 

amongst  the  various  molecular  species  present  to  possess  a  permanent  electric 

dipole  moment,  it  seems  reasonable  to  attempt  an  explanation  of  the  tensile 
stress  in  terms  of  the  mutual  repulsion  of  molecular  dipoles  adsorbed  on  the 
walls  of  pores  known  to  exist  within  the  film.  On  this  basis  one  would  expect 
the  generation  of  a  tensile  stress  not  to  be  a  property  specific  to  MgF2  only. 
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but  one  that  is  shared  by  other  substances  with  an  appropriate  porous  structure, 
in  particular  by  materials  such  as  cryolite  and  ZnS  that,  in  common  with  MgF2 
have  columnar  habits  of  growth.  From  the  work  of  Pulker  and  Jung (ref. 5), 
Kinosita  and  Nishibori(ref.6)  and  Ogura  et  al(ref.7)  it  is  known  that  these 
substances  have  pore  structures  that  are  readily  penetrated  by  atmospheric 
moisture.  The  development  of  tensile  stress  has  not  yet  been  investigated  in 
these  other  materials,  but  our  work  has  shown  that  thin  films  of  shellac* 
develop  tensile  stress  in  a  similar  manner,  and  of  a  magnitude  comparable  to 
that  found  in  MgF2 . 

In  contrast  to  this,  vacuum  deposited  silver  films  do  not  exhibit  the  effect 
at  all,  although  they  are  readily  permeated  by  water  vapour.  This  permeability 
can  be  demonstrated  by  depositing  Ag  onto  a  MgF2  film  and  observing  the 
deflection  of  the  composite  system  on  exposure  to  water  vapour.  Figure  3  shows 
the  results  of  such  an  experiment,  in  which  a  500  A  silver  layer  was  put  down  on 
the  same  MgF2  film  that  had  been  the  subject  of  the  earlier  experiments  relating 
to  figure  2,  and  since  preliminary  measurements  has  established  that  the  elastic 
properties  of  the  arrangement  were  not  significantly  affected  by  the  additional 
silver  deposit,  direct  comparison  of  the  data  of  figures  2  and  3  is  possible. 

In  forming  the  Ag  film,  care  was  taken  to  ensure  that  no  portion  of  the  MgF2  was 
accessible  to  H20  other  than  by  penetration  of  the  silver  barrier,  and  bearing 
in  mind  that  the  latter  on  its  own  did  not  develop  any  tensile  stress,  the 
subsequent  movement  of  the  cantilever  must  be  ascribed  to  water  vapour  that  had 
passed  through  the  barrier.  In  comparing  this  beam  deflection  with  that  shown 
in  figure  2,  it  is  evident  that  whilst  the  rate  of  pressure  rise  is  similar  in 
both  cases,  the  high  impedance  offered  by  the  silver  film  to  the  passage  of 
H20  leads  to  a  far  more  sluggish  response  of  the  compound  film,  so  that  the 
equilibrium  beam  deflection  is  reached  only  about  two  hours  after  the  steady 
state  pressure  level  has  been  attained.  It  will  also  be  noticed  that,  although 
this  pressure  level  is  higher  than  in  the  earlier  experiment  of  figure  2,  it 
results  in  a  smaller  ultimate  deflection  of  the  beam.  We  believe  this  may  be 
due  to  some  of  the  pores  in  the  MgF2  being  blocked  off  by  the  silver  film,  and 
thus  rendered  inaccessible  to  the  H2O  vapour. 

The  fact  that  any  water  vapour  at  all  could  penetrate  to  the  MgF2  clearly 
implies  a  degree  of  porosity  in  the  silver  layer,  but  to  account  for  the  absence 
of  any  detectable  tensile  stress  when  the  latter  on  its  own  is  exposed  to 
moisture,  we  suggest  that  the  number  of  pores  per  unit  area  of  the  Ag  film  may 
be  too  low.  Alternatively  the  average  pore  diameter  could  be  too  small  to 
accommodate  the  number  of  polar  molecules  required  to  produce  a  significant 
tensile  force.  This  we  shall  discuss  more  fully  in  the  next  section,  where  we 
develop  a  simple  model  to  describe  how  the  tensile  stress  is  related  to  the 
number  of  molecular  dipoles  adsorbed  in  the  porous  structure  of  the  film. 


3.  THEORY  OF  THE  ADSORPTION  STRESS 

In  constructing  a  model  for  the  stress  mechanism  we  are  not  at  this  stage 
concerned  with  developing  a  theory  of  great  numerical  precision.  Our  aim  is 
rather  to  establish  whether  or  not  an  explanation  in  terms  of  dipole  inter¬ 
actions  between  adsorbed  molecules  is  at  all  feasible,  i.e.  whether,  on  the 
basis  of  such  interactions,  the  observed  tensile  forces  can  be  produced  by 
realistic  amounts  of  the  adsorbate.  In  our  particular  case,  involving  water 
vapour  at  a  pressure  of  a  few  Torr,  adsorbed  on  a  surface  at  ambient  temperature, 
we  take  this  to  imply  a  surface  coverage  roughly  in  the  range  from  one  to  ten 
monolayers. 

Following  Pulker  and  Jung  (loc.cit.)  we  represent  the  columnar  growth 
structure  of  the  MgF2  film  by  an  aggregate  of  closely  packed  cylindrical  columns 
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Prepared  by  painting  mica  substrates  with  a  solution  of  1  g  of  shellac 
flakes  in  15  cc  of  methyl  alcohol. 
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of  radius  R  (see  figure  4(a)).  In  this  model  "pores"  of  cross-sectional  area 
(\/3  -  tt/2)R2  are  formed  between  any  three  adjacent  cylinders.  To  simplify  the 
calculations,  we  replace  the  resulting  rather  complex  pore  shape  by  circularly 
cylindrical  pores  of  radius  a,  chosen  so  as  to  yield  the  same  cross-sectional 
area. 

From  simple  geometrical  considerations  it  is  found  that  these  "equivalent" 
pores  have  a  radius: 


0.23  R 


(2) 


and  that  the  spacing  between  adjacent  pore  centres: 


(V§~  -  1A^3  a 


(3) 


We  now  consider  the  adsorbed  molecular  dipoles  to  be  arranged  on  the 
cylindrical  pore  wall  with  their  axes  arranged  normal  to  the  surface,  and  with 
charges  of  the  same  sign  pointing  inward  and  outward  respectively.  We  shall 
take  the  circumference  of  the  adsorbing  surface  to  remain  at  the  constant  value 
of  27Ta,  irrespective  of  the  amount  of  material  adsorbed,  i.e.  we  disregard  the 
decrease  in  area  available  for  adsorption  as  the  thickness  of  the  adsorbate 
layer  becomes  appreciable. 

Let  the  length  of  the  pore  be  equal  to  the  film  thickness  T^,  and  let  us 

introduce  a  cartesian  co-ordinate  system,  the  z-axis  of  which  is  coincident  with 
the  pore  axis,  whilst  the  x-and  y-axes  lie  in  the  centre  plane  of  the  film,  with 
the  x-axis  pointing  along  the  length  of  the  mica  strip  i.e.  in  the  direction  of 
the  tensile  stress  (see  figure  4(b)).  A  dipole  of  moment  Mj  is  aligned  along 
the  x-axis  in  the  x-y  plane,  and  we  wish  to  examine  its  interaction  with  a 
second  dipole  of  moment  M2 ,  arranged  at  a  vertical  height  h  above  the  centre 
plane  and  displaced  azimuthally  from  Mi  by  an  angle  <j>.  The  force  along  vector 
r  connecting  Mi  and  M2  is  given  from  elementary  electrostatic  theory(ref .8)  as: 


^r  =  ^efr|  4  (sin#  sin#1  cos^  -2COS0COS#1)  (4) 

Here  e  is  the  dielectric  constant  of  free  space,  ^  is  the  angle  between  the 
planes  intersecting  in  r  that  contain  the  dipole  vectors  Mi  and  M2 ,  whilst  # 
and  01  are  the  angles  between  r  and  the  two  dipole  axes,  as  indicated  in  the 
figure. 

Introducing  the  normalised  vertical  height: 


y 


(5) 


it  can  readily  be  shown  that: 


sin# 


sin2p+  y2 


(2sin(0/2))2  +  sin20  +  y2 


(6) 
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sin#1  = 


sin2  0(4  sin2  (0/2)  (sin2  (0/2)  +  cos0)  +  1)  -  sin40  +  y2 
(2sin4  (0/2))2  +  sin20  +  y* 


(7) 


COS#  = 


2 sin2  (0/2) 


V(2sin2  (02))2  +  sin20  +  y2 


(8) 


cos 


i  _  2sin2  (0/2)cos0-sin20 

V(2sin2  (0/2) )2  +  sin20  +  y2 


(9) 


and 


COS0  = 


y2  cos 0  +  sin 20 
y2  +  sin 20 


(10) 


In  addition  the  absolute  value  of  the  vector  r  is  given  by: 

I  rl  =  a  V4sin2  (0/2)  +  y2  (11) 

Equations  (6)  to  (11)  express  the  various  quantities  appearing  in 
equation  (4)  as  functions  of  0  and  of  the  normalised  vertical  height  y. 
Therefore,  making  the  appropriate  substitutions,  the  interaction  force  acting 
along  r  is  obtained  in  terms  of  these  variables  only.  Multiplication  by 
F 

cos  0  yields  x  Mi  ,  the  x-  component  of  the  interaction  force  on  Mi  ,  which 
represents  the  contribution  to  the  tensile  force  we  are  concerned  with,  and 
which  from  equations  (4)  and  (6)  to  (11)  may  be  written  in  the  form: 


xFMi 


5MiMz 
47re  a4 


g(0>y) 


(12) 


where  g  is  a  function  of  0  and  y  only. 

We  require  the  tensile  force  on  Mi  due  not  only  to  its  interaction  with  M2 , 
but  due  to  that  with  all  the  polar  molecules  adsorbed  on  the  pore  wall. 

Rather  than  carrying  out  a  summation  over  the  entire  set  of  discrete  molecular 
dipoles,  it  is  more  convenient  for  the  present  to  replace  these  by  an 
equivalent  continuous  dipole  layer  of  moment  P  per  unit  area,  and  to  perform  an 
integration,  replacing  Mi  and  M2  by  dipole  elements  P  dSi  and  P  dS2  of  area  dSi 
and  dS2  respectively.  In  our  approximate  calculation  this  procedure  is 
permissible  since,  as  shown  by  Topping (ref. 9)  the  potential  at  a  point  in  a 
dipole  array  is  rather  insensitive  to  a  rearrangement  of  the  individual  dipoles, 
provided  the  change  is  made  in  a  manner  that  leaves  the  average  dipole  moment 
per  unit  area  constant. 

Noting  that: 


dS2  =  a2  d0dy 


the  total  component  of  the  tensile  force  on  dMi  then  becomes: 
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3/3 2  dS! 
Aire  a2 


,,  .  ...  3/32dSi  T 

g(0,y )  d^dy  =  - —  I 

47rea2 


(14) 


The  value  of  the  interaction  integral  I  obviously  depends  on  the  integration 
limits  of  the  variables  4>  and  y,  and  it  is  these  limits  that  need  further 
consideration.  The  angular  variable  </>  clearly  has  the  upper  limit  it ,  and  for 
a  dipole  element  dMi  on  the  midplane  of  the  film  y  ranges  up  to  T^/2,  but  for  an 

element  not  situated  on  this  plane  we  assign  a,  for  the  moment,  unspecified 

upper  limit  y  .  Also,  as  a  first  step  we  take  the  lower  limits  <j>  and  y  to 
J  max  o  o 

be  both  zero.  Values  of  the  integral  I  computed  on  this  basis  are  shown  in 

the  curve  of  figure  5  as  a  function  of  y  .  It  is  seen  that  at  y  =  0.5  I 
°  max  max  o 

has  practically  reached  its  asymptotic  value  of  about  8.66.  Since  h  =  ay,  this 

means  that  dipoles  separated  from  dMi  by  a  vertical  distance  h  in  excess  of  half 

a  pore  radius  make  only  a  small  contribution  to  the  interaction  integral.  In 

practice  T^  ^  a,  and  we  can  therefore  always  set  ymax  =  °°  >  irrespective  of 

whether  or  not  dMi  lies  in  the  midplane  of  the  film.  In  doing  so  we  introduce 
only  an  edge  effect  error,  involving  the  small  fraction  of  molecules  adsorbed 
on  the  pore  wall  within  a  distance  of  less  than  a/2  from  the  film  boundary. 

By  setting  the  lower  limits  and  yQ  equal  to  zero  we  are  allowing  the 

continuous  dipole  distribution  to  approach  the  axis  of  the  dipole  element  dMi 
indefinitely  closely.  This  leads  to  an  overestimate  of  the  interaction  force, 
as  in  actual  fact  two  neighbouring  molecular  dipoles  are  separated  by  a  minimum 
distance  d  ,  representing  in  our  case  the  effective  diameter  of  the  water 
molecule.  °  To  obtain  a  better  approximation  we  therefore  restrict  the  closest 
approach  to  dQ  by  setting  =  dQ/a  for  0  <  y  <  dQ  and  =  0  for  y  >  dQ. 

Equation  (14)  then  assumes  the  form; 


/Ml  = 

27rea2 


f 

J 


7T 


do/a 


g(tf,y)#dy  + 


7T 


g(0,y)d$dy 


5/3 2  dSi 
2ttc  a2 


I (a,  dQ) 


(15) 


where  a  factor  2  has  been  introduced  to  take  into  account  the  contribution  from 
points  both  above  and  below  the  centre  plane. 

For  a  dipole  element  dMi  not  oriented  along  the  x-axis,  but  inclined  to  it 
at  an  angle  <£,  the  contribution  to  the  tensile  force  acting  on  the  film  cross- 
section  is  obtained  by  multiplying  equation  (15)  by  cosV.  For  the  entire  pore 
wall  one  thus  finds: 


3/32  I  (a,dQ) 


-7T/2 


6/32  TfI  (a,dQ) 


F 

x  pore 


ite  a 


co  s'P  d^dh  = 


7rea 


(16) 
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Since  there  are  1/s  —  l/10a  pores  per  unit  length  (see  equation  (3)),  the 
total  tensile  force  acting  on  the  cross-section  of  the  film  per  unit  width  is 
given  by: 


total 


602I(a,dQ)Tf 

lCwrea2 


(17) 


Equating  this  to  the  total  tensile  force  as  derived  from  the  observed  beam 
deflection  by  means  of  equation  (1)  and  solving  for  0  one  obtains  for  the 
dipole  moment  per  unit  area: 


Da  /  20  7r  e  E  5 
L  V  9i(a,do)T£ 


(18) 


At  this  stage  it  is  convenient  to  abandon  the  model  of  an  equivalent 
continuous  dipole  layer,  that  we  introduced  earlier  for  ease  of  calculation, 
and  instead  regard  jS  more  realistically  as  being  due  to  n  dipoles  per  unit 
area,  each  having  a  dipole  moment,  a,  so  that: 

0  =  nd  (19) 


and  the  number  density  of  polar  molecules  adsorbed  on  the  pore  surface  is 
obtained  as: 


n 


Da  J  20jreE5__ 
«L  V  9  I(a,do)Tf 


(20) 


In  order  to  evaluate  this  expression  we  must  introduce  appropriate  numerical 
values  for  the  effective  molecular  diameter  dQ  and  for  the  pore  radius  a. 

From  the  work  of  Pulker  and  Jung  (loc.cit.)  we  expect  the  pore  radius  to  lie 
roughly  in  the  range  from  10  A  to  20  A,  and  the  molecular  diameter  may  be 
estimated  from  the  area  A  covered  by  a  molecule  on  the  substrate,  an  area  that 
we  shall  take  to  be  circular.  Following  Dushman(ref . 10)  we  write  this  area  as: 


A 


4  \/3 


M 


,4  \/2  N  p 

ci 


2/3 


(21) 


where  N  is  Avogadro's  number,  whilst  p  and  M  are  the  density  and  molecular 

cL 

weight  of  the  condensed  adsorbate.  The  effective  molecular  diameter  is  then 
found  as: 


d  =  1.4  x  10“8  (M/  )1/3  cm  (22) 

yielding  for  H20  a  diameter  of  3.6  A.  This  represents  an  appropriate  value 
only,  but,  as  seen  in  figure  6(a),  changes  of 'the  order  of  10%  to  20%  in  dQ  do 

not  affect  the  interaction  integral  I(a,d  )  in  a  manner  that  would  invalidate 

our  analysis.  Also,  according  to  figure  6(b),  for  dQ  =  3.6  A 
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I(a,dQ)  s  5.8  x  108a 


(23) 


where  a  is  expressed  in  metres. 

Equation  (20)  can  therefore  be  written  as: 


n  —  3.3  x  10“10  ^jr  -\j  molecules/m2 


(24) 


Introducing  for  our  particular  case: 


beam  thickness  D  = 
beam  length  L  = 
Young's  Modulus  E  = 
film  thickness  = 

dipole  moment  a  = 
pore  radius  a  = 

we  obtain: 

n  = 


1.6  x  10”5  m 
2  x  10“2  m 
1.8  x  1010  kg/m2 
4  x  10”7  m 

6.1  x  10'30  M.K.S.  units 
1  x  10-9  m 

2.9  x  1020  \/T  molecules/m2 


(25) 


For  5  =  5  x  10-4  m  this  yields  n  —  6.5  x  1018  molecules/m2 ,  representing 
roughly  one  monolayer,  a  value  that  probably  underestimates  n  by  a  factor  of 
order  unity,  an  underestimate  that  arises  at  least  in  part  from  the  fact  that 
we  have  neglected  any  decrease  in  the  dipole  moment  u  due  to  mutual  depolarisa¬ 
tion  of  neighbouring  molecules.  In  addition  we  have  also  selected  a  pore 
radius  in  the  lower  part  of  the  expected  range.  Nevertheless  the  calculated 
value  for  n  is  of  a  reasonable  order  of  magnitude (ref .11,12) ,  thus  confirming 
that  an  explanation  of  the  tensile  stress  in  terms  of  dipole  interactions  is 
feasible. 

Equation  (25)  establishes  a  relation  between  n  and  the  observed  beam 
deflection  5.  On  the  other  hand,  from  general  adsorption  theoryfref . 13)  there 
is  a  further  relation  connecting  n  with  the  adsorption  energy  Q  through  the 
equation: 


n  =  3’5  p  r  eQ/RoT  molecules/m2  (26) 

Vmt 

where  Rq  is  the  gas  constant  (cal/mole) ,  T  the  absolute  temperature,  and  r  is 

the  vibration  period  of  the  adsorbed  particles,  to  which  we  assign  the  value 
of  1  x  10~13  s  (de  Boer,  loc.cit.  p35) .  In  the  above  conventional  form  of  the 
equation  Q  is  regarded  as  a  constant,  an  approach  that  is  permissible  when  the 
adsorption  takes  place  on  a  homogenous  surface,  and  the  surface  coverage  is 
sufficiently  low  to  exclude  interaction  between  the  adsorbate  molecules,  but  in 
our  case  it  is  more  appropriate  to  replace  Q  by  Q(n) ,  a  suitably  weighted 
average  value  over  all  adsorbed  molecules,  and  one  that  will  in  general  vary 
with  n.  By  equating  expressions  (25)  and  (26)  we  then  find  for  T  =  300  K  that: 


v'l  .  1.65  x  10-’ a^"^600 

P 


(27) 
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The  ratio  v  5/p  can  be  derived  from  our  experimental  data  (see  figure  2),  and 
as  the  graph  of  figure  7  shows,  it  may  be  satisfactorily  represented  over  the 
range  of  our  measurements  by  the  empirical  equation: 

VJ  _  9.9  x  IQ"13 

iT  "  5 3  <28) 

where  5  is  expressed  in  metres  and  p  in  Torr.  Substituting  in  equation  (27) 
we  obtain  the  average  adsorption  energy  in  terms  of  the  beam  deflection  as: 


Q(n) 


600  l  n 


(29) 


The  solid  portion  of  the  curve  in  figure  8  represents  the  average  adsorption 
energy  calculated  in  this  manner,  but  it  must  be  realised  that,  because  of  its 
reliance  on  the  empirical  formula  (equation  (28)),  equation  (29)  has  only  a 
limited  range  of  applicability,  particularly  since  our  measurements  of  the 
dependence  of  5  on  p  extend  only  up  to  8  s  6  x  10“4  m,  and  for  beam  deflections 
below  about  2  x  10-4  m  the  data  are  only  qualitative  due  to  the  limited 
accuracy  of  the  manometer  in  the  corresponding  pressure  range. _  This  accounts 
for  the  fact  that  equation  (29)  predicts  an  unlimited  rise  of  Q(n)  as  5  0, 

when  in  fact  we  must  expect  it  to  approach  the  adsorption  energy  for  the  bare 
adsorbent  surface  in  the  manner  indicated  qualitatively  by  the  upper  dotted 
portion  of  the  curve.  Similarly  a;t  large  values  of  5,  when  a  number  of 
adsorbate  layers  have  been  formed,  Q(n)  should  gradually  decrease  to  a  value 
near  the  heat  of  vaporisation  of  water  in  the  condensed  phase  (  1 0  000  cal/ 
mole) .  This  asymptotic  value  is  seen  to  be  quite  compatible  with  the  calculated 
curve  of  figure  8,  although,  in  view  of  the  various  numerical  uncertainties  in 
our  calculation,  such  a  high  degree  of  compatibility  must  to  some  extent  be 
fortuitous.  Nevertheless  the  fact  that  the  calculated  adsorption  energies  lie 
in  the  range  expected  from  earlier  experiments  pn  water  vapour  adsorption 
(de  Boer,  loc.cit.  p42) ,  lends  further  credibility  to  our  model. 

As  the  total  surface  area  of  the  pore  system  within  the  film  is  covered 

by  adsorbate  molecules  to  a  surface  density  n  (in  our  example  6.5  x  1018 
molecules/m2)  an  amount  of  energy: 


nA  Q(n) 

H  =  p  (30) 

N 

a 


is  released,  and  the  question  arises  as  to  whether  this  energy  is  sufficient 
to  account  for  the  energy  of  elastic  deformation  stored  in  the  cantilever.  A 
simple  calculation  shows  that  under  our  conditions  this  energy  is  about 
5  x  10-9  cal  for  a  beam  deflection  of  5  x  10-4  nu_  The  same  value  of  5  yields 
from  equation  (29)  an  average  adsorption  energy  Q(n)  of  10  500  cal/mole.  Also, 
with  the  aid  of  equation  (3)  the  area  of  the  pore  system  in  a  film  of  width  W 
(in  our  case  4  x  10-3  m)  is  given  by: 


2TTfWL  x  10'2 


(31) 


so  that  for  a  =  1  x  10-9  m  one  obtains  that  H  =  2.3  x  10-4  cal.  The  energy 
that  becomes  available  during  the  adsorption  process  therefore  exceeds  the 
stored  elastic  energy  by  more  than  four  orders  of  magnitude. 
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4.  CONCLUSION 

The  results  of  the  previous  sections  show  that  the  tensile  forces  produced  in 
a  thin  porous  film  by  the  adsorption  of  water  vapour  can  be  accounted  for  in 
terms  of  the  electrostatic  interaction  of  molecular  dipoles,  and  according  to  our 
model  polar  molecules  other  than  HjO  should  cause  similar  effects.  Moreover  it 
is  apparent  from  equations  (17)  and  (23)  that  the  forces  due  to  the  dipole  inter¬ 
action  are  approximately  inversely  proportional  to  the  pore  size,  so  that  little 
stress  will  be  found  in  films  with  very  coarse  porosity.  The  converse  prediction 
of  high  stress  when  the  pore  size  is  small  does  however  need  some  qualification, 
since  our  analysis  assumes  the  thickness  of  the  adsorbate  to  be  small  in 
comparison  with  the  pore  radius.  This  assumption  will  not  in  general  hold  for 
very  fine  pores,  where  moreover  the  number  of  adsorbed  molecules  will  not  only  be 
determined  by  the  dynamic  equilibrium  of  adsorption  (equation  (26))  but  an  upper 
limit  to  this  number  will  in  addition  be  set  by  the  amount  of  adsorbate  that  can 
be  accommodated  within  the  confines  of  the  small  pore  volume.  As  a  result,  both 
very  small  and  very  large  pores  will  lead  to  insignificant  adsorption  stress, 
but  there  is  an  intermediate  range  of  pore  radii,  placed  by  our  analysis  roughly 
between  10  A  and  100  A,  where  very  substantial  stresses  (in  our  experiments  of 
the  order  of  2  x  107  kg/m2 )  can  be  induced. 

If  this  porosity  range  can  be  avoided  by  a  suitable  choice  of  film  deposition 
parameters,  there  should  not  only  be  a  considerable  reduction  in  adsorption 
stress,  but  it  would  seem  that  also  the  optical  stability  of  the  films  could  be 
markedly  improved.  It  has  long  been  recognised  that  the  refractive  index  n  of 
MgF2  films  changes  following  exposure  to  atmospheric  moisture.  In  the  past  this 
effect  has  been  ascribed  to  water  filling  the  voids  within  the  film,  and  in  doing 
so  raising  the  refractive  index  of  the  pore  space  from  1.0  for  air  to  1.33  for 
H20,  thus  increasing  the  average  refractive  index  of  the  film  by  an  amount 
depending  on  the  total  pore  volume.  On  the  basis  of  this  picture  several  authors 
have  in  fact  deduced  the  film  porosity  from  changes  in  the  value  of  n. 

On  the  other  hand  when  a  stress  S  is  applied  to  a  material,  its  refractive 
index  can  also  be  changed  through  the  photo-elastic  effect  by  an  amount: 

A  n  =  BS  (32) 

where  B  is  the  appropriate  stress-optical  coefficient,  measured  in  Brewsters,  if 
S  is  expressed  in  dynes/cm2 (ref . 14) .  It  is  therefore  of  interest  to  examine 
whether  the  stress  induced  by  H2O  adsorption  could  materially  contribute  by  this 
mechanism  to  the  observed  variation  in  the  refractive  index. 

A  stress  of  about  2  x  109  dyn/cm2 ,  as  encountered  in  the  present  work,  together 
with  a  typical  refractive  index  change  A  n  of  5  x  10  3 (ref. 6),  would  imply  for  B 
a  value  of  25  Brewsters.  Whilst  no  photo-elastic  data  for  MgF2  appear  to  be 
available,  B  =  25  falls  into  a  typical  range  for  a  number  of  materials (ref. 14,15) , 
thus  making  our  hypothesis  at  least  feasible. 

We  can  obtain  further  supporting  evidence  from  the  manner  in  which  n  changes 
with  the  moisture  content  of  the  atmosphere.  Since  from  equation  (1)  S  is 
linearly  related  to  the  beam  deflection  5,  whilst  according  to  equation  (28)  § 

is  in  turn  proportional  to  p2^7,  the  photo-elastic  change  in  n  (equation  (32)) 
can  be  expressed  in  terms  of  the  partial  pressure  p  of  the  water  vapour  as: 

A  n  =  BlP2/7  (33) 

A  variation  in  accordance  with  this  power  law  is  in  fairly  good  agreement 
with  experimental  observations,  as  shown  in  figure  9 (a), -where  data  are 
reproduced,  obtained  by  Kinosita  and  Nishibori(ref .6)  for  a  MgF2  film  exposed 
to  an  atmosphere  of  variable  moisture  content,  the  latter  being  expressed  as  a 
fraction  of  the  saturated  vapour  pressure  pQ.  The  change  in  refractive  index 
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is  seen  to  be  associated  with  hysteresis,  and  in  the  figure  the  solid  curves  are 
as  drawn  by  the  authors,  whilst  the  dotted  lines  were  calculated  from  an 
expression  of  the  form  of  equation  (33)  with  Bj  set  equal  to  7.3  x  10~3  .  This 
value  was  arrived  at  by  fitting  the  power  law  to  the  lower  curve  of  figure  9(a), 
although  acceptable  agreement  with  the  upper  curve  is  also  obtained,  in  spite  of 
the  fact  that  the  latter,  in  view  of  the  small  number  of  experimental  points 
given,  is  only  semi-quantitative  in  nature. 

The  agreement  of  equation  (33)  with  the  experimental  data  of  Kinosita  and 
Nishibori  provides  further  support  for  our  suggestion  that  adsorption  induced 
stress  may,  via  the  photo-elastic  effect,  contribute  substantially  to  the 
observed  changes  in  the  refractive  index.  Should  this  be  the  case,  earlier 
estimates  of  film  porosity  based  on  variations  in  n(ref.6,7)  may  need 
reassessment. 

It  will  however  be  noticed  that  significant  departures  from  the  2/7  power  law 
occur  for  values  of  p/PQ  in  excess  of  about  0.8.  It  is  believed  that  these 

deviations  may  be  due  to  the  effect  of  capillary  adsorption  (de  Boer,  loc.cit. 
p210) ,  and  such  an  interpretation  could  also  explain  a  similar  departure  in  the 
same  range  of  p/pQ  (see  figure  9(b)),  found  by  Kinosita  and  Nishibori  in  the 

case  of  a  MgF2  film  exposed  not  to  H20  but  to  ethyl  alcohol  vapour,  a  substance 
with  a  dipole  moment  comparable  to  that  of  water.  The  figure  also  shows  that 
for  p/p  below  0.8  the  measurements  are  very  convincingly  represented  by  a  2/7 
power  law,  suggesting  that  our  empirical  equation  (28) ,  from  which  the  power  law 
is  derived,  is  of  greater  generality  than  might  initially  have  been  expected. 
Thus  the  equation  is  seen  to  yield  a  good  approximation  for  observations  made  on 
a  number  of  films,  prepared  in  different  laboratories  and  exposed  to  different 
polar  adsorbents.  Particularly  in  view  of  this  latter  fact,  it  seems  that  the 
equation  reflects  a_  rather  general  trend  in  the  manner  in  which  the  average 
adsorption  energy  Q(n)  changes  as  the  number  of  adsorbed  molecules  increases. 
Because  of  its  fundamental  implications  this  aspect  will  need  further 
investigation. 

Turning  to  some  more  immediately  practical  applications,  the  compound  film 
technique  (Ag/MgF2 )  discussed  in  Section  2  affords  a  ready  means  of  determining 
the  penetration  of  H2O,  and  presumably  of  other  polar  molecules,  through  thin 
films,  provided  the  pore  structure  of  these  films  is  such  as  to  prevent 
significant  adsorption  stress  in  the  material,  the  permeability  of  which  is 
under  study. 

On  the  other  hand  films  of  materials  such  as  MgF2 ,  that  exhibit  large 
adsorption  stresses,  constitute,  in  combination  with  a  suitable  elastic  sub¬ 
strate,  the  sensing  element  of  a  simple  and  very  sensitive  hygrometer.  Initial 
experiments  point  to  a  detection  limit  of  such  a  device  in  the  vicinity  of 
5  x  10-4  Torr  partial  pressure  of  H20  (~  3  x  10-3%  relative  humidity),  with  good 
reproducibility  up  to  pressures  of  several  Torr,  although  there  is  some 
hysteresis  when  the  elements  are  first  exposed  to  very  high  water  vapour 
pressures.  This  appears  to  be  eliminated  however  by  repeated  cycling  through 
the  humidity  range.  The  long  term  performance  of  such  elements  is  now  being 
investigated,  and  will  be  reported  in  due  course. 
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Figure  2.  Deflection 


BEAM  POSITION  (mm) 


Figure  3.  Effect  of  water  vapour  on  a  MgF2  film  covered  by  500  A  of  Ag 
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Figure  7.  \/8/p  as  function  of  1/53 


Figure  8.  Mean  adsorption  energy  as  function  of  beam  deflection 


Figure  9.  Effect  of  polar  molecule  adsorption  on  the  refractive  index  of  MgF: 
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